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We report the experimental and the analytical results of electron spin resonance (ESR) 
and magnetization in high magnetic fields up to about 68 T of the quasi two-dimensional 
triangular lattice antiferromagnet NiGa2S4. From the temperature evolution of the ESR 
absorption linewidth, we find a distinct disturbing of the development of the spin correlation 
by ^2-vortices between 23 K and 8.5 K. Below T v = 8.5 K, spin-wave calculations based on 
a 57° spiral spin order well explains the frequency dependence of the ESR resonance fields 
and high field magnetization processes for H\\c and 7J_Lc, although the magnetization for 
H-Lc at high fields is different from the calculated one. Furthermore, we explain the field 
independent specific heat with Independence by the same spin-wave calculation, but the 
magnitude of the specific heat is much less than the observed one. Accordingly, these results 
suggest the occurrence of a Zi vortex-induced topological transition at T v and may indicate 
quantum effects beyond the descriptions based on the above classical spin models. 

KEYWORDS: high field ESR, high field magnetization, triangular lattice, frustrate magnetism, 
NiGa 2 S 4 



1. Introduction 

Two-dimensional (2D) triangular lattice antiferromagnets (TLAFM) have attracted con- 
siderable attention over the three decades, because some of them are believed to have a novel 
ground state such as a spin liquid caused by strong frustration. 1 ) Recently, a new quasi-2D 
S=l TLAFM NiGa2S4 has attracted a substantial interest as a candidate of such a spin liq- 
uid state. 2 - 1 It has slabs consisting of two GaS layers and a MS2 one separated by a van der 
Waals gap as shown in Fig. 1(a). Despite strong antiferromagnetic (AF) interactions with 
Weiss temperature 9 W ~— 80 K, no long-range order (LRO) down to 0.08 K was observed. 3 ) 
Some noteworthy characteristic features were observed in neutron diffraction, specific heat, 
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and magnetic susceptibility measurements. 2 ) The neutron diffraction measurements unveil the 
development of a short-range correlations with the wavevector q = (0.158, 0.158, 0), which 
correspond to a short-range order (SRO) with about 57° spiral spin structure as shown in Fig. 
1(b). The correlation length of the SRO is only about 2.5 nm (about seven lattice spacings) 
even at 1.5 K, and the origin of the absence of LRO has not been clarified yet. The magnetic 
susceptibility (M/H) measured at 0.01 T exhibits a weak peak at about 8.5 K and a finite 
value at T — > K, indicating the existence of gapless excitations. The magnetic specific heat 
shows broad peaks at about 10 K and 80 K, and T 2 -dependence with no field dependence up 
to 7 T below 10 K. From the T 2 -dependent specific heat, one may associate the existence of 
a A;-linear excitation in a 2D system attributed to a Nambu-Goldstone mode resulting from a 
breaking of a continuous symmetry. These curious experimental results have stimulated further 
theoretical and experimental studies on this compound. 3 ~ 14 * > A theoretical analysis by taking 
bilinear-biquadratic exchange interactions into consideration seemed to account for some of 
such peculiar experimental results. 4 ~ 6 ) Strong third-neighbor and negligible second-neighbor 
exchange interactions are suggested by analyses of a photoemission spectroscopy 7 ' and ab 
initio calculations by a density functional theory (DFT). 8 ) Taking the strong third-neighbor 
exchange interactions into account, a finite-temperature phase transition is discussed using a 
Monte Carlo method. 9 ) Recent NMR and //SR measurements evidenced the occurrence of a 
novel state where the spins have a MHz-scale dynamics through a weak freezing below about 
10 K, 3 ' 11 ) and impurity effect studies of Nii_ a; M a ;Ga2S4 (M = Mn, Fe, Co and Zn) suggested 
that the integer size of the Heisenberg spins is important to stabilize the 2D coherent behavior 
in the frozen spin-disordered state. 12 ) 

In this paper, we report the experimental and analytical results of electron spin resonance 
(ESR) and magnetization measurements on single crystals of NiGa2S4, which were performed 
to clarify the spin dynamics in more details than those in our previous short report. 13 ) We 
found peculiar dynamics of Z<i vortex-induced effect on the temperature dependence of ESR 
absorption linewidth. The frequency dependence of the ESR resonance fields, high field mag- 
netizations and low-temperature specific heat below 8.5 K are explained to a some extent 
by a conventional spin-wave theory, which is consistent with spin-wave excitations expected 
in the Zi vortex-induced phase. These results suggest an occurrence of a Z2 vortex-induced 
topological transition at 8.5 K. However, some results are quantitatively different from the 
calculated ones based on a classical spin model. In particular, the calculated specific heat 
value below about 10 K is largely different from the observed one, suggesting an additional 
quantum effects such as quadrupolar correlations. 
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2. Materials and Methods 

2.1 Crystal structure 

First we summarize the structural properties of NiGa2S4. Single crystals were grown by 
chemical- vapour transport using iodine. The details of crystal growing was reported in Ref. 
14. NiGa2S4 belongs to a trigonal system with P3ml symmetry and the lattice constants a = 
b = 0.3624 nm and c = 1.1999 nm, and has a layered unit consisting of central magnetic layer 
of edge-sharing NiS6 octahedra and nonmagnetic two layers of GaS4 tetrahedra, in which the 
magnetic layer is sandwiched, as shown in Fig. 1(a). These layered units are separated from 
each other along the c-axis by a van der Waals gap, resulting in highly 2D system with a weak 
ferromagnetic interlayer interaction. 2 ) In the N1S2 layers, Ni 2+ (3d 8 ) ions with an electronic 
configuration t\ g e 2 g (S = 1) form the exact TLAFM as shown in Fig. 1(b). 

2.2 Experimental methods 

ESR measurements in pulsed magnetic fields up to about 53 T were carried out at temper- 
atures between 1.3 K and about 100 K using our pulsed field ESR apparatus equipped with a 
non-destructive pulse magnet, a far-infrared laser (Edinburgh Instruments, U.K.) which cover 
the frequencies up to about 1.4 THz. We also performed ESR measurements at 1.4 K in static 
magnetic fields up to 14 T utilizing a superconducting magnet (Oxford Instruments, U.K.) 
and a vector network analyzer with some extension (ABmm, France) for the frequencies up to 
about 500 GHz. High-field magnetization measurements in pulsed magnetic fields up to about 
68 T were carried out with a non-destructive pulse magnet at 1.3 K at ISSP in The University 
of Tokyo. The magnetization was measured by an induction method using a pick-up coil. In 
the experiments, external fields were applied parallel (H _L c) and perpendicular (H || c) to 
the 2D planes. 

3. Experimental Results 

3.1 High field ESR 

Observed ESR absorption spectra are displayed in figures 2(a)-(d). Figures 2(a) and 2(b) 
show the temperature dependence of the ESR absorption spectra at 584.8 GHz for H \\ c 
and H _L c, respectively. We observed a strong resonance signal, whose linewidth is extremely 
broad, and a weak resonance signal for both directions. The resonance fields are indicated 
by the arrows. Since the weak resonance signals are almost temperature independent, it is 
considered to arise from paramagnetic resonances caused probably by imperfections in the 
lattices. The temperature dependences of the resonance field and the absorption linewidth of 
the strong signal are shown in Figs. 3 and 4, respectively. The resonance field shifts to a lower 
field with decreasing temperature for both directions below about 30 K. The inset of Fig. 3 
shows the g- values converted from the resonance fields. The absorption linewidth increases as 
the temperature is lowered from about 80 K, corresponding to the absolute value of the Weiss 
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temperature \6 W \. The temperature dependence of the linewidth shows a bending at about 23 
K, below which the increasing tendency with decreasing temperature becomes gradual. Below 
about 8.5 K, the linewidth shows a gradual decrease as the temperature is lowered. 

The frequency dependences of the ESR absorption spectra at 1.3 K for H || c and H 
_L c are presented in Figs. 2(c) and 2(d), respectively. The resonance fields are indicated by 
the arrows. All the resonance fields are plotted in the frequency- field plane as seen in Fig. 
5. The resonance modes from the strong signals have a zero-field gap and come close to the 
paramagnetic resonance line in the higher field for both H || c and H _L c, and such behavior 
of the ESR modes is similar to that observed in some spiral antiferromagnets with easy- 
plane anisotropy. 15 ) Since the resonance frequencies of the weak resonance signals are almost 
proportional to the external fields, they can be considered to arise from the paramagnetic 
resonance as mentioned in the previous paragraph. 

3.2 High field magnetization 

Figures 6(a) and 6(b) indicate magnetization curves at 1.3 K in pulsed magnetic fields 
for H || c and H _L c, respectively. For the former, the magnetization increases linearly with 
increasing field for the whole range of applied magnetic field, while, for the latter, it deviates 
from a linear increase shown by the broken line above about 30 T. In addition, we observed 
a small anomaly at about 44 T as shown in the inset of Fig 6(b). 

4. Analyses and Discussion 

4-1 Temperature dependence of the ESR linewidth 

From the temperature dependence of the ESR linewidth as shown in Fig. 4, we divide the 
temperature region into three areas, (i) above 23 K, (ii) between 23 K and 8.5 K, and (iii) 
below 8.5 K. In the high-temperature region (i), we analyze the temperature dependence of 
the ESR linewidth in terms of (a) a common expression of the 2D critical behavior and (b) 
a non-linear sigma (NLcr) model 16 ) associated with the AF LRO in a 2D Heisenberg AFM. 
In the case (a), the temperature dependence of the linewidth is expressed by the following 
equation, 

AH 1/2 oc (T — r N )"f , (1) 

where Tn is the AF long-range ordering temperature, p is a critical exponent, which reflects 
the anisotropy and the dimensionality of the system. Based on Kawasaki's dynamic equations, 
the critical exponent of the 2D AFM is given as p = (3 — 2rj)v } where v is the critical exponent 
for the inverse of the spin correlation length, n the critical exponent related to the static spin 
correlation which is assumed to obey the form, <S*Sq> = (k 2 + q 2 )^ l+r il 2 } 7 ) In the case 
of the 2D Heisenberg AFM, p is evaluated to be about 2.5 from the experimentally obtained 
T] = 0.2 and v = 0.95, and the value p is known to be close to the experimental values in 
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several Heisenberg 2D compounds. 18 -* In NiGa2S4, we expect that the development of the SRO 
gives the same critical behavior as in such 2D Heisenberg AFMs. Since there is no LRO down 
to 0.08 K, 3 ) we assume that Tn is approximate to zero. Then, we fit Eq. (1) with p = 2.5 and 
Tn — to the experimental linewidth and obtain good agreement between them as shown by 
the dotted line in Fig. 4. 

In the case (b), we use the following relation of the temperature dependence of the ESR 
linewidth calculated by Chakravarty and Orbach for an S = 1/2 nearest-neighbor quantum 
Heisenberg AFM on a square lattice with a small interplanar coupling, 

+ (2) 

where £ is a spin-correlation length, k-Q is the Boltzmann constant, and p s is a spin-stiffness 

constant. 16 ^ The temperature dependence of the spin-correlation length for 2D frustrated 

quantum antiferromagnets predicted by a NL<r model is given by 19 ) 

1 ,4vrps , . 

£°c^^exp(- — — ). (3) 

The p s /k-B, is evaluated to be about 18 K from a NLcr model analysis of dynamic relaxation 
rates of /xSR measurements, u ) which is the same value as that evaluated from the specific 
heat. 2 ) We fit Eq. (2) to the experimental results with this value. The theoretical result in 
Fig.4 (gray line) indicates a more rapid increase of the linewidth than the experimental one. 
This disagreement may originate from strong spin frustration as discussed in the study of the 
nuclear spin-lattice relaxation rate in Ref. 20. 

Below 23 K, the increase of the linewidth with decreasing temperature becomes gradual, 
indicating that some effect disturbs the development of the SRO below 23 K. In the temper- 
ature region (ii), we examine two types of stable point defect in 2D systems, (c) a Z2 vortex 
for the Heisenberg TLAFM and (d) a conventional Kosteritz-Thouless (KT) vortex for XY 
TLAFM. The Z2 vortex is a topologically stable point defect characterized by a two-valued 
topological quantum number. 21 ) Unlike the KT vortex of spins, it is a vortex of spin chirali- 
ties. The vortices form bound pairs at low temperatures and begin to dissociate at a certain 
critical temperature Ty like a KT-type phase transition. Above Ty, the unbound free vortices 
are thermally excited, and the spin fluctuates through the spacial passage of the vortices. The 
effect of the Z2 vortices on the ESR linewidth was discussed on Heisenberg TLAFMs, HCr02 
and LiCr02- 22 ) According to the analyses in these compounds, the ESR linewidth above Ty 
is expected to be written as 

Ai7 1/2 ocexp(g;), (4) 

where Ey is the activation energy of a free Z2 vortex. A Monte Carlo simulation gives Ty = 
0.31J5 2 and Ey = 1.65JS 12 , 21 ) and thus the ratio Ey /Ty becomes constant. We take Ty = 8.5 
K from the anomaly of the magnetic susceptibility and the peak of the ESR linewidth and 
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then obtain Ey = 45.2 K. Using this value, satisfactory agreement between experiment and 
calculation is obtained between 8.5 and 23 K as shown by the solid line in Fig. 4. In the neutron 
diffraction experiment, it was found that the correlation length in the triangular plane grows 
up to be about three lattice units at 23 K. 23 ) Our result indicates that a correlation length 
of about three lattice units is needed to observe the Z<i vortex effect. The shift of the ESR 
resonance field below about 30 K is also thought to be caused by the development of the 
correlation length which is about two lattice units at 30 K. 23 ) 

For the case (d), it is expected that thermally excited free KT vortices fluctuate the spin 
system above a certain critical temperature Tkt- Neglecting the temperature dependence of 
the vortex velocity, the temperature dependence of the ESR linewidth was reported to be 
proportional to a cube of the correlation length and given by 24 ) 

37T 

AH 1/2 xf = exp( ). (5) 

The fit obtained with Tkt = 8.5 K cannot reproduce the gradual slope of the experimental 
result as shown by the dotted gray line in Fig. 4. The KT-type phase transition is accompanied 
by a divergence of the correlation length at Tkt- 25 ' 1 in NiGa2S4, no divergence was observed in 
such experiments. This is consistent with the fact that the obtained temperature dependence 
of the ESR linewidth deviates from Eq. (5). 

In the low-temperature region (iii), the ESR linewidth shows a gradual decrease as the 
temperature is lowered and it remains a large value even at the lowest temperature. If there is 
a conventional phase transition to a LRO phase at about 8.5 K, the ESR linewidth is expected 
to decrease rapidly with decreasing temperature reflecting a development of uniform internal 
fields. 26 ) In the case of NiGa2S4, from the NQR measurements, Takeya et al. suggested that a 
gradual spin freezing occurs between 10 K and 2 K and that a inhomogeneous internal fields 
appear below 2 K. 3 ) A considerably large ESR linewidth at the lowest temperature must 
reflect a distribution of the inhomogeneous internal fields, which is consistent with the NQR 
result. 



4-2 ESR resonance modes in the low-temperature region 

Next, we discuss the frequency dependence of the ESR resonance fields of the strong signals 
at the lowest temperature. Figures 5(a) and 5(b) show the resonance fields at various frequen- 
cies for H || c and H _L c, respectively. The broken lines represent paramagnetic-resonance 
lines with each g factor, = 2.08 and g± = 2.17. The ESR modes of the strong signal have a 
zero-field energy gap and approach the paramagnetic-resonance modes with increasing field for 
both field directions. These modes are reminiscent of the spiral spin resonance mode observed 
in a TLAFM with easy-plane anisotropy such as RbFe(Mo04)2- 15 ^ In addition, it is known 
that some low-dimensional magnets in the temperature region where the SRO well develops 
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exhibit resonance modes expected in the LRO phase. 27 ) Therefore, we analyze the experimen- 
tal results in terms of conventional spiral spin resonances. Considering strong third-neighbor 
and negligible second-neighbor exchange interactions, the calculation based on a mean-field 
approximation shows that the exchange constants are required to be J1/J3 = —0.2 and J 2 = 
to realize the 57° spiral spin structure indicated by the neutron diffraction measurements. 
The Ji (i = 1, 2, 3) represent the first-, second-, and third-neighbor exchange interactions, 
respectively. Then, the spin Hamiltonian is written as, 

n = Ji Si-Sj + Sl ' Sm + D Y,^ 2 ~ 9» B J2 s i- H ( 6 ) 

<ij> <lm> i i 

where D is the anisotropy of the easy-plane type (D > 0), fiB the Bohr magneton, H the 
external magnetic field, the z-axis parallel to the c-axis, (ij) all the nearest neighbor pairs, 
and (Im) all the third-neighbor pairs. In our previous study, for simplicity, we assumed J\=0, 
which results in a 120° spiral spin structure on the third nearest neighbor pairs in the 
plane. 13 ) In the present study, we take into account the finite J\ to realize the 57° spiral spin 
structure, which is formed by 120 sublatices. Then, the free energy F is expressed by the 
following form using the mean-field approximation, 

„ 120 120 

F = A Mi-Mj + M v M m + -^(M? ) 2 - ^M, H, (7) 

<ij> <lm> i=l i=l 

where A, B, and K is given by, 

,120 Ji R= 120_7^ 120 2D 

N (gps)*' N ( 5/UB ) 2 ' N ( gf i B )^ W 

and Mi is the i-th sublattice moment and given by 

N 

Mi = —g^ B S l , (9) 

where TV is the number of magnetic ions, and Si is the spin on the z-th sublattice. We derive 
the resonance conditions by solving the equation of motion 

dMi/dt = 1 [M i xH i }, (10) 

where 7 is the gyromagnetic ratio and Hi a mean-field applied on the i-th sublattice moment 
given by 

Hi = -dF/8M t . (11) 

To solve the equation of motion, we use a method applied for ABX3 type antiferromagnets. 28 ) 
Assuming precession motions of the sublattice moments around those equilibrium directions, 
we utilize the following expressions, which represent the motion of the i-th. sublattice moment, 

Mi = (AM i:i exp(iujt),AM^exp(iujt),\Mi\). (12) 
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where AM^, AM^<^.\Mi\, and x, y and z are the principal axes of the coordinate system 
on each sublattice moment. The i-axis is defined to be parallel to the direction of the each 
sublattice moment, and the x and y-axes are perpendicular to that. 

In the case of H || c, all the sublattices are tilted from the easy plane with equivalent angle 
as shown by the schematic picture in Fig. 6(a). The angle 9 between each individual sublattice 
moment and the external field is determined from the equilibrium condition dF/dO = 0, and 
9 is given by 



where (j) = 57°. Considering the change of 9 under the external field, we obtain 120 resonance 
modes by solving Eq. (10) numerically. If time dependent parts, namely transverse compo- 
nents, of the sublattice moments cancel each other and thus the total magnetization is not 
in motion, the resonance modes are, in principle, not observable. From the calculation of the 
total magnetization, we obtain two observable gapped modes, hv\ and hu2, and a gapless 
mode, /if 3 = 0. The hu\ calculated with the parameters Ji/k-Q = —4.56 K, J^/k-Q = 22.8 K, 
D/k-Q = 0.8 K exhibits good agreement between the experiment and the calculation as shown 
in Fig. 5(a). The zero-field energy gap is evaluated to be about 13 K, which indicates a sub- 
stantial anisotropy and decreases with increasing temperature to be about 10 K at T v = 8.5 
K. Since the hv-i has a weak field dependence, it is expected to show a much broader resonance 
linewidth than hv\, thus being hardly observed. In the case of H _L c, we assume the distorted 
spin structure, in which spins lie in the easy plane as shown by the schematic view in Fig. 
6(b). The distortion angles ipi are determined so as to minimize the energy of the assumed 
structure by solving differential equations dF/dcpi = (i = 1~120), and we acquire them as 
numerical solutions. Then, through the same procedure as in the case of H || c, we obtain 
two gapped modes, hv 4 and hi/5, and a gapless mode, hv 6 = 0. The hi/4 mode shows good 
agreement with the experimental one using the same parameters for H \\ c as shown in Fig. 
5(b). However, we did not observe the signals for hv§. The reason is that the signal intensity 
of hv§ mode is expected to be very weak due to the canceling of the transverse components of 
the sublattice moments, which is crucial for ESR transitions. This mode becomes soft around 
45 T where one third of the saturation magnetization occurs. 

Here, we briefly discuss the realization of a Z2 vortex-induced topological transition on 
our system with a considerably large D value. Quite recently, Kawamura derived the relation 
D/kB<Ty(^/a)~ 2 which is a necessary condition for the Z2 vortex-induced transition. 29 ) 
In the case of a classical TLAFM model with only nearest-neighbor interaction, the spin 
correlation length £ is estimated to be O(10 3 ) at T = T v so that even an extremely small D 
would destroy the Z2 vortex transition. However, if the spin correlation length £ is supressed 
due to some reason, the anisotropy is of the same order of T v (£/a) -2 or less, and the Z2 



cos 9 = 



g\\VBH 



(13) 



2S(D + 3(Ji + J 3 ) - 2 Ji cos <f> - ( Ji + 2 J 3 ) cos 2<j) - J 3 cos 40) ' 
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4-3 Magnetization curves in the low-temperature region 

Using the parameters determined by the ESR analyses, we calculated the magnetization 
curves. For H _L c, we take into account not only a spiral structure but also a so called 
umbrella structure with an out-of-plane spiral spin configuration. Then, from a mean-field 
theory, a first-order phase transition from the spiral structure to the umbrella structure is 
expected to occur at about 25 T. No such behavior, however, is observed and thus, spins must 
lie in the easy plane even above 25 T for H _L c. The magnetizations are written as 

My =g llt i B S cos d, (14) 

120 

M ± = g ± v B S^^. (15) 
i=i 

As shown in Fig. 6(a), the agreement between the experiment and the calculation is excellent in 
magnetic fields up to 63 T for H \\ c. We reproduce the non- linear increase of magnetization 
for H _L c as shown in Fig. 6(b) by considering continuous changes of ipi. The calculation 
suggests a first-order phase transition from a spiral spin structure to a fan structure at about 
63 T for H _L c, where the spins oscillate sinusoidally along the external field direction. 
However, the first order transition was not observed in our experiment. We cannot explain 
this disagreement within the framework of the mean-field theory. In a Heisenberg TLAFM with 
the nearest-neighbor AF interactions, a phase transition from a spiral spin configuration to a 
collinear spin one, in which the spins on the two sublattices are parallel and those on the other 
sublattice are antiparallel to the applied magnetic field in the basal plane, is expected to occur 
at one third of the saturation magnetization due to quantum or thermal fluctuations. 30 ' 31 - ) 
A small anomaly of the magnetization observed at 44 T where the magnetization is close 
to one third of the saturation value, 0.72 /ie/ion, mav be related to the 1/3 magnetization 
plateau discussed in the TLAFM. In our mean- field approximation corresponding to a classical 
calculation at T = 0, the collinear configuration exists at a singular field, 31 ) and the calculated 
magnetization does not show any 1/3 magnetization plateau. 

4-4 Specific heat in the low-temperature region 

The above ESR resonance modes obtained by the mean-field theory correspond to the 
spin-wave excitations as a function of magnetic field at q = 0, ±0.158, where q = ak/2-ir (k is 
a wavevector and a is the lattice constant of the triangular plane). This suggests that spin- 
wave-like excitations develop in the basal plane at low temperatures despite the absence of the 
LRO in NiGa2S4. In Fig. 7, we show spin-wave dispersion energies, hu k , at H=0 T and H=7 
T for H || c calculated by means of a conventional spin- wave theory with parameters obtained 
by our ESR analyses. We find a fc-linear dispersion around k ~ 0, which remains gapless even 



9/19 



J. Phys. Soc. Jpn. Full Paper 

in a magnetic field up to 7 T. Furthermore, the spin-wave velocity v s , which corresponds to 
the slope of the /c-linear dispersion, is almost constant in magnetic fields up to 7 T. Actually, 
we obtain v s = 1045 m/s at 7 T and v s = 1049 m/s at T, which are almost equivalent to each 
other and are nearly the same value as that evaluated from a neutron scattering experiment. 23 ) 
This almost field-independent dispersion must cause T 2 -dependent specific heat below about 
10 K. Therefore, we evaluate the specific heat in the low-temperature region. The specific heat 
is derived from the derivative of the average thermal energy written as 

C = dE/dT. (16) 

From the classical spin-wave theory, the E is given by 

E = ^ D(k)dk, (17) 

J o exp(^) - 1 

where D(k) is a density of states. In the present triangular lattice model, it is given by, 

D(k) = ^±k. (18) 

where Na is the Avogadro number. Below about 10 K, we only consider the linear dispersion 
around k ~ because the energy gap was evaluated to be about 13 K, and the spin-wave 
dispersion obeys the relation hv^v s k. Accordingly, the specific heat obtained from the above 
relations is calculated as 

C=^N A kU^f, (19) 

where £(3) = 1.202, h Planck's constant. Although we confirm that the obtained v s is close 
to the value evaluated from the recent neutron diffraction measurement, 23 ) the calculated C 
value is nearly 36 times smaller than the experimental one. For quantitative agreement with 
the experimental specific heat, the value of v s is required to be nearly six times 32 -* smaller 
than that obtained from our ESR analyses. We consider that the quantitative disagreement 
in the present spin-wave model may come from some kind of effect, for example, anomalous 
structure factor of the excitation mode or quadrupolar correlations stabilized by a biquadratic 
interaction. 4 ' 5 ^ For the latter, Kawamura and Yamamoto predicted by means of their Monte 
Carlo simulations that any type of biquadratic interaction could couse a Z2 vortex-induced 
topological transition. 6 ^ 

5. Conclusions 



We performed high field ESR and magnetization measurements on single crystals of the 
triangular lattice antiferromagnet NiGa2S4. We find three distinct temperature regions: (i) 
above 23 K, (ii) between 23 K and 8.5 K, and (iii) below 8.5 K. In the high-temperature 
region (i), the temperature dependence of the ESR linewidth indicates a development of the 
SRO of the 2D Heisenberg AFM. In the intermediate temperature region (ii), the development 
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of the SRO is suppressed by thermally excited Z2 vortices and the temperature dependence of 
the ESR linewidth is proportional to the exponential function with the activation energy of the 
free Zi vortex toward 8.5 K. In the low-temperature region (hi), the frequency dependence of 
the ESR resonance fields and the magnetization processes are explained to a large extent by the 
conventional mean-field theory. Furthermore, we derive a gapless linear dispersion relation for 
H || c from the conventional spin- wave theory and explain qualitatively the field- independent 
specific heat with T 2 dependence. The origin of the large quantitative disagreement between 
experiment and calculation based on classical spin models suggest significant contribution 
from quantum effects, possibly, due to quadrupolar correlations. Although there are some 
distinctions from the conventional spin-wave model, the observed properties, explained by the 
spin- wave excitations in spite of the absence of LRO, are consistent with those expected in 
the Z2 vortex- induced phase. 6 ) Consequently, unconventional properties of NiGa2S4 in the 
low-temperature region strongly suggest that a vortex-induced topological transition takes 
place at 8.5 K. 
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Fig. 1. (Color online) (a)Structure of the unit slab of NiGa2S4 consisting of two GaS layers and one 
MS2 layer stacking along the c-axis. (b)Short range correlated spiral spin structure of NiGa 2 S4 at 
low temperature region. The a&-plane components of spins rotate with an angle about 57°. 
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Fig. 2. (Color online) (a), (b) Temperature dependence and (c), (d) frequency dependence of ESR 
absorption spectra of NiGa 2 S 4 at 584.8 GHz and at 1.3 K for H \\ c ((a),(c)) and H _L c ((b), (d)). 
The closed and open arrows indicate the broad large signal and the weak one, respectively. The 
sharp signals at higher field come from an ESR standard sample of DPPH for correction of the 
magnetic field. 
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Fig. 3. (Color online) Temperature dependence of the resonance field of the strong ESR signals in 
NiGa2S4 at 584.8 GHz. Open and closed circles represent the resonance fields for H \\ c and H _L 
c, respectively. The inset shows g- values converted from the resonance fields. 
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(b) NLa model 
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Fig. 4. (Color online) Temperature dependence of the full width at half maximum of the broad ESR 
spectra in NiGa2S4 for H \\ c at 584.8 GHz. Black broken, gray solid, black solid, and gray broken 
lines represent the temperature evolution of the linewidth calculated for (a) 2D critical behavior, 
(b) non-linear sigma (NLer) model, (c) ^-vortex scenario, and (d) Kosterlitz-Thouless (KT) vortex 
scenario, respectively. The inset shows the expansion of the temperature region between 8 K and 
23 K. 
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Fig. 5. (Color online) Frequency-field plot of the resonance fields of NiGa2S4 measured at 1.3 K in 
pulsed magnetic fields and at 1.4 K in static magnetic fields for (a) H \\ c and (b) H _L c. Closed 
circles and crosses denote strong and weak signals, respectively. Solid and broken lines represent 
calculated resonance modes and the paramagnetic-resonance lines, respectively. 
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Fig. 6. (Color online) Magnetization curves of NiGa2S4 at 1.3 K for (a) H || c and (b) H _L c. Solid 
and blue broken lines are experimental and calculated magnetization curves, respectively. The 
black thin broken line is an extension of the liner increase at low-field region. The inset of (b) 
shows the vicinity of the anomaly for flic. The illustration in each figure shows a schematic 
view of spin configuration under an external magnetic field for each direction. 
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Fig. 7. (Color online) Calculated spin-wave dispersion relation of NiGa2S4. Solid and broken lines 
represent the relations at H = T and H = 7 T for H || c, respectively. 
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